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We report on the quasi-two-dimensional magnetism of the natural mineral diaboleite 
Pb2Cu(OH)4Cl2 with a tetragonal crystal structure, which is closely related to that of the frustrated 
spin- 1 magnet PbVC>3. Magnetic susceptibility of diaboleite is well described by a Heisenberg spin 
model on a diluted square lattice with the nearest-neighbor exchange of J ~ 35 K and about 5 % 
of non-magnetic impurities. The dilution of the spin lattice reflects the formation of Cu vacancies 
that are tolerated by the crystal structure of diaboleite. The weak coupling between the magnetic 
planes triggers the long-range antiferromagnetic order below Tjv — 11 K. No evidence of magnetic 
frustration is found. We also analyze the signatures of the long-range order in heat-capacity data, 
and discuss the capability of identifying magnetic transitions with heat-capacity measurements. 
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I. INTRODUCTION 

Quantum spin systems show intricate low-te mper ature 
phenomena of fundamental^ and even applied^ inter- 
est. Despite strong quantum fluctuations that tend to 
impair and eventually destroy ordered spin states, most 
quantum magnets develop long-range magnetic order at 
sufficiently low temperatures. In two and three dimen- 
sions, the lack of the ordered state at zero temperature 
is only possible in frustrated magnets, where the com- 
petition between magnetic couplings amplifies quantum 
fluctuations. For example, systems of interest are those 
based on spin-| in the square-lattice or kagome-lattice 
geometriesP The simple square lattice with nearest- 
neighbor Heisenberg interactions (Ji) is non- frustrated. 
However, a second-neighbor coupling Ji frustrates the 
system and leads to a rich phase diagram of the J\ — J 2 
frustrated square lattice model that was extensively stud- 
ied in the pastPHm 

Experimental search for the J\ — J 2 systems on a 
square lattice remains a challenging problem. While 
materials with ferromagnetic (FM) J\ can be prepared 
in a rather systematic fashion using the building blocks 
of V +4 phosphates] 12 * 13 ^ systems with antiferromagnetic 
(AFM) Ji are less studied. In the limit of J2/J1 — > 0, 
Cu +2 -based coordination compounds^^ and VOM0O4 
(Ref. IT8|) are excellent material prototypes of square- 
lattice systems with weak magnetic frustration. The 
stronger frustration with J2/J1 approaching 0.5 has been 
so far observed only in one material, PbVC>3 (J2/J1 — 
0.35), that remains controversial because muon spin ro- 
tation experiments detect the magnetic ordering below 
43 Kp^O but neither thermodynamic measurements nor 
neutron powder diffraction are capable of observing this 
magnetic transition! 2 ^ Latest theoretical results suggest 
that above J%jJ\ = 0.35 — 0.4 the Ji — J 2 square-lattice 
system enters the spin-liquid regime, and the magnetic 
order vanishes While PbVC>3 may lie on (or be close 
to) the boundary of this spin-liquid region, further sys- 
tems with J2/J1 — 0.5 are highly desirable. 



The crystal structure of PbVC>3 is a tetragonal deriva- 
tive of the perovskite type. The tetragonal distor- 
tion, which is essential for the quasi-two-dimensional 
(2D) square-lattice magnetism, is only rarely observed 
in perovskites. Similar structures are found in BiCoC>3 
(Ref. |2"T|) and in the "supertetragonal" polymorph of 
BiFe03, which is stabilized in thin films. 22 However, none 
of these compounds is suitable as a model system for the 
quantum spin-| square-lattice model because both Co +3 
and Fe +3 bear much higher and largely classical spins. 
Surprisingly, an appropriate structural analog of PbVC>3 
can be found in a completely different family of materi- 
als. The Cu-based mineral diaboleite, Pb 2 Cu(OH)4Cl2, 
has tetragonal crystal structur e! 23 * 24 ! that can be derived 
from the perovskite structure type. 

The structure of diaboleite features isolated CUO4 pla- 
quettes that form magnetic layers in the ab plane. These 
layers are interleaved by Pb2Cl2 slabs (Fig. [I]). The re- 
lation to the parent ABX3 perovskite structure can be 
understood by writing the chemical composition of dia- 
boleite as Pb 2 (CuD)X 6 , where X = OH, CI. This way, 
half of the B-type positions are vacant, with the result- 
ing voids filled by the bulky OH~ anions. While this 
non-trivial transformation leaves little chemical similar- 
ity to perovskites, the crystal structure of diaboleite still 
features the space group Pimm, same as in PbV0 3 and 
PbTiC>3. This tetragonal crystal structure is perfectly 
suited for a quasi-2D magnetism, because the magnetic 
layers of Cu +2 atoms are well separated, with the short- 
est interlayer distance of 5.5 A. The robust tetragonal 
symmetry entails the perfect square-lattice arrangement 
of the magnetic sites in the ab plane. 

In the following, we present an experimental and com- 
putational study of diaboleite. We show that this com- 
pound is indeed a good material prototype of the spin- g 
square lattice. However, its J%j J\ ratio is very small, 
hence the magnetic frustration is weak. We discuss the 
origin of this low J2, and pinpoint another peculiarity, 
the dilution of the spin lattice through the formation of 
Cu vacancies that may have implications for experimen- 
tal studies of the diluted square lattices with spin- ~ . 
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FIG. 1. (Color online) An overall view of the diaboleite crystal structure (left panel) and its projection on the ab plane (right 
panel). Large spheres depict the Pb (violet) and CI (brown) atoms located between the [Cu(OH)4] magnetic layers. Green 
spheres inside the plaquettes are Cu atoms. Small spheres show the O (orange) and H (gray) atoms. The middle panel shows 
the natural sample of diaboleite used in the present study. Note the translucent blue crystals of diaboleite and opaque white 
crystals of laurionite and phosgenite (see Sec. Ill] for details). 



The outline of this paper is as follows. After summariz- 
ing experimental and computational methods in Sec. [TTJ 
we report the detailed characterization of the natural 



mineral sample (Sec. Ill A I and elaborate on possible de- 



viations from the ideal crystal structure (Sec. MB). We 
further develop microscopic magnetic models for both the 
ideal and distorted crystal structures (Sec. |III C ), and 
apply them to thermodynamic properties of diaboleite 
(Sec. HID). Our results are discussed in Sec. IV which 



is followed by a summary and an outlook. 



II. METHODS 

A natural sample of diaboleite (Mammoth- Saint 
Antony Mine, Pinal Co., Arizona, USA) was provided 
by the mineralogical collection of the Department of Ma- 
terials Research and Physics at Salzburg University. A 
visual inspection of this sample (Fig. [I] middle panel) 
identified translucent blue crystals of diaboleite mixed 
with opaque white crystals of other minerals. Both blue 
and white crystals were detached from the support and 
separated manually. This way, a 100 mg batch of the 
diaboleite crystals was obtained. 

Phase composition of the samples was determined by 
laboratory x-ray diffraction (Huber G670 Guinier cam- 
era, CuK a i radiation, ImagePlate detector, 29 = 3 — 100° 
angle range). The batch of the blue crystals mostly 
contained diaboleite, whereas white crystals were iden- 
tified as a mixture of phosgenite (Pb2Cl2C03) and lau- 
rionite (PbOHCl). Considering possible contamination 
of the natural samples with variable - both crystalline 
and amorphous - impurities, sample quality was further 
checked by bulk chemical analysis and high-resolution 
x-ray diffraction (XRD). The chemical analysis was per- 
formed with the ICP-OES method.^ The high-resolution 
XRD data were collected at the ID31 beamline of the Eu- 
ropean Synchrotron Radiation Facility (ESRF, Grenoble) 
at the wavelength of about 0.43 A. Details of the experi- 
ment are described elsewhere.^ The Jana2006 program 
was used for the structure refinement P3 Crystal struc- 



tures were visualized using the VESTA software.^ 

Magnetization measurements were performed with 
Quantum Design MPMS SQUID magnetometer in the 
temperature range 2 — 380 K in fields up to 5 T. Heat 
capacity was measured with Quantum Design PPMS in 
the temperature range 1.8 — 100 K in zero field and in 
the fields of 3 T, 6 T, and 9 T. All measurements were 
performed on polycrystalline samples. 

Magnetic couplings in diaboleite, as well as differ- 
ent aspects of its structural distortions, were analyzed 
by density- functional theory (DFT) band-structure cal- 
culations performed in the FPL(PH and VASP^l codes 
that implement basis sets of local orbitals and pro- 
jected augmented waves respectively. Local density 
approximation (LDA^P^and generalized gradient approx- 
imation (GGA)32l for the exchange-correlation potential 
were used. Band-structure results are obtained for well- 
converged k meshes with 336 points in the symmetry- 
irreducible part of the first Brillouin zone for the crys- 
tallographic unit cell, and about 100 points for the su- 
percells. Residual forces in optimized crystal structures 
were below 0.01 eV/A. Effects of strong electronic cor- 
relations in the Cu 3d shell were treated by introducing 
the LDA band structure into an effective Hubbard model 
(model approach), or within the mean-field DFT+£7 
procedure, with the on-site Coulomb repulsion param- 
eter Ud = 6.5 eV and the Hund's exchange parame- 
ter Jd — 1 eV applying the around-mean-field double- 
counting correction schemePS Details of the computa- 
tional procedure are described in Sec. |III C| 

Thermodynamic properties and Neel temperatures of 
the spin models relevant to diaboleite were obtained from 
quantum Monte-Carlo (QMC) simulations based on the 



loop^and dirloop_sse' :15 algorithms of the ALPS^ sim- 
ulation package. Simulations were performed on 2D 
(L x L) and 3D (Lxix L/5) finite lattices with L < 80 
and L < 50, respectively. Size convergence for the ther- 
modynamic properties (magnetization and magnetic spe- 
cific heat) was carefully checked. The magnetic ordering 
temperatures (TV) were determined by calculating tem- 
perature dependence of the spin stiffness p s for different 
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FIG. 2. (Color online) High-resolution XRD pattern of the 
natural sample of diaboleite. Ticks show the reflection posi- 
tions, whereas asterisks mark tiny reflections of the phosgenite 
impurity. The reflections of the second impurity, wherryite, 
are not visible in this angular range. Note that the 003 reflec- 
tion is narrower than the reflections with non-zero h and/or k. 



L and applying the appropriate scaling procedure, as pre- 
viously described in Refs. [37] and [35] 



III. RESULTS 
A. Sample characterization 

The high-resolution XRD pattern of the natu- 
ral diaboleite sample reveals strong reflections of 
the tetragonal diaboleite phase along with few very 
weak reflections that are attributed to the im- 
purities of phosgenite (Pb2Ci2C03) and wherryite 
[Pb 7 Cu2(S0 4 )4(Si04)2(OH) 2 ]. The amounts of these im- 
purities are 1.1(1) wt.% and 0.4(1) wt.%, respectively, ac- 
cording to the Rietveld refinement.^ While phosgenite is 
a diamagnet, no information on the magnetism of wher- 
ryite is presently available. Nevertheless, trace amounts 
of this impurity should not affect any of the results pre- 
sented belowPS 

Chemical analysis delivered the following bulk com- 
position of the natural sample: 68.1(1) wt.% Pb, 
9.3(1) wt.% Cu, and 0.2(1) wt.% Ba. Other elements 
amenable to the ICE-OES determination (including tran- 
sition metals, but excluding O, H, and CI) are below 
the detection limit of 0.1 wt.%. These estimates are in 
reasonable agreement with the expected chemical com- 
position (67.2 wt.% Pb, 10.3 wt.% Cu), although slight 
deviations in the Pb and Cu content, as well as the pres- 
ence of Ba, signal subtle variations in the stoichiometry 
of the diaboleite phase. 

The Rietveld refinement of the atomic positions re- 
sulted in reasonable atomic displacement parameters 
(ADPs) of Pb, CI, and O, U iso ~ 0.01 A 2 (the hydro- 
gen position was not refined) However, the notably 
higher ADP of Cu (U iso ~ 0.035 A 2 ) suggested the miss- 



ing scattering density at the Cu position. Considering 
the nearly isotropic thermal ellipsoid (U^ n ~ 0.03 A 2 , 
U^? u ~ 0.04 A 2 along the a and c directions, respectively), 
the large ADP of Cu can not be explained by a displace- 
ment of Cu from its position on the four-fold axis, and 
should be rather understood as the formation of vacan- 
cies. Indeed, as the ADP was fixed to Ui SO = 0.01 A 2 , the 
refined Cu occupancy converged to 0.942(4) indicating 
the deficiency of Cu atoms in the diaboleite structure.^ 
The formation of Cu vacancies conforms to the underesti- 
mated Cu content from the chemical analysis that yields 
about 10 % Cu deficiency in the diaboleite sample. 

The formation of Cu vacancies is justified by both pow- 
der XRD and chemical analysis. Indeed, the missing scat- 
tering density in XRD could be otherwise explained by 
the substitution of a lighter element in the Cu position. 
However, the chemical analysis does not show detectable 
amounts of any foreign elements that are capable of re- 
placing Cu in the diaboleite structure. Trace amounts of 
Ba are likely related to the Pb/Ba substitution, which is 
favored by the similar ionic radii of these elements. The 
amount of the substituted Ba is well below 1 %, so it 
does not show up in the XRD refinement. 

Another conspicuous feature of the natural diaboleite 
sample is the sizable broadening of the hkl reflections 
with non-zero indices h and/or k. This effect is well seen 
in Fig. [2] where, for example, the 310 reflection is much 
broader than the neighboring 003 peakpO Anisotropic re- 
flection broadening can be understood as the formation 
of stacking faults in the layered crystal structure of dia- 
boleite. This problem is further addressed in Sec. |III B[ 
where we also elaborate on the nature of Cu vacancies in 
the diaboleite structure. 



B. Details of the crystal structure 

The structure refinement puts forward several issues 
regarding details of the atomic arrangement that may 
be important for understanding the magnetism of dia- 
boleite: i) positions of hydrogen atoms that are not pre- 
cisely determined by XRD; ii) formation of Cu vacancies; 
iii) stacking faults that underlie the reflection broaden- 
ing. 

The position of hydrogen has been determined from 
single-crystal XREB31 using the soft constraint on the 
O-H bonds (the O-H distance of 0.98 A). We further 
refined the hydrogen position by a DFT-based struc- 
ture optimization, with all atoms other than hydrogen 
fixed to their experimental positions. Both LDA and 
GGA optimizations^ converged to very similar results. 
The hydrogen atom is in the 4<i position (x, x, z) with 
x = 0.2975/0.2972 and z = 0.1191/0.1179 in LDA/GGA, 
respectively. The resulting O-H distance is 0.99 A, and 
the Cu-O-H angle is about 112.3°. These parameters 
only slightly differ from those found experimentally: the 
O-H distance of 0.98 A (constrained) and the Cu-O-H 
angle of 109.1°. 
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FIG. 3. (Color online) Left panel: formation of Cu vacancies 
in the magnetic layer of diaboleite (4 of Cu atoms and the 
respective amount of oxygens are removed). Top right panel: 
the pristine structure of stoichiometric diaboleite. The re- 
moval of one Cu and one O atom releases a water molecule. 
Additionally, the coordination of Pb atoms changes from four 
longer Pb-O bonds (2.45 A, top right panel) to two shorter 
Pb*-0 bonds (2.32 A and 2.40 A, respectively; left panel). 
Only the Pb-O bonds shorter than 2.5 A are shown. The 
notation of atoms follows Fig. [T] 



The formation of Cu vacancies can be explained as 
follows. The removal of one Cu atom creates an uncom- 
pensated negative charge that should be balanced by the 
removal of one oxygen or two chlorines. As the oxy- 
gen atoms are strongly bonded to Cu with the Cu-0 
distance of 1.97 A, the removal of oxygen is expected. 
The remaining hydrogen atom may join one of the OH 
groups to form a water molecule. Regarding the four 
oxygen atoms that surround the vacancy, one is removed 
to compensate the charge, one oxygen forms the water 
molecule, whereas the two remaining OH groups are still 
bound to Pb. The resulting structure has been fully re- 
laxed within GGA+U to yield the configuration shown 
in Fig. [3] where we consider the 2x2 supercell with one 
out of four Cu atoms removed. This supercell with 25 % 
Cu deficiency should be representative for the actual dia- 
boleite sample that features 5—10 % of Cu vacancies. The 
calculations for lower concentrations of vacancies would 
require even larger supercells, which are hardly feasible 
for an accurate computational treatment. 

The formation of Cu vacancies enhances the Pb-O 
bonding and reduces the coordination number of Pb 
atoms around the vacancy. In stoichiometric diaboleite, 
each Pb atom is surrounded by four oxygens with the 
Pb-O separations of 2.45 A. In Cu-deficient diaboleite, 
the Pb* atoms (Fig. [3]) are strongly bound to only two 
O atoms, but the respective Pb-O distances are short- 
ened to 2.32 A and 2.40 A, respectively. The bond to the 
third oxygen atom extends to 2.53 A, whereas the fourth 
bond is either lost (for Pb* atoms adjacent to the oxy- 
gen vacancy) or stretched up to 2.63 A (for the oxygen 
atom of the water molecule). This way, the lengths of 
the Pb-O bonds are redistributed, and bonding require- 
ments of Pb can be satisfied, even though one out of four 
oxygen atoms is removed from the crystal structure. The 




FIG. 4. (Color online) Possible mechanism of stacking 
faults in the diaboleite structure. The notation of atoms fol- 
lows Fig. [l] Left panels: parent tetragonal crystal structure. 
Right panels: optimized structure with the uniform shift of 
the [Cu(OH) 4 ] layers for k = (||0). Note that the [Pb 2 Cl 2 ] 
slabs are invariant to this translation, which is denoted by 
the thick arrow. The weak CI • • H hydrogen bonds are the 
only interaction between the [Pt>2Cl2] slab and the neighbor- 
ing [Cu(OH)4] layer. Therefore, the structures with k = 
(left panel) and k = (||0) (right panel) marginally differ in 
energy. In the bottom panels, hydrogen atoms are omitted 
for clarity. 

void formed between the two Pb* atoms (Fig.[3| is typical 
for Pb +2 oxides. It can be ascribed to the "localization" 
of 6s 2 lone pairs that form a non-bonding region in the 
crystal structure.^ 

The formation of Cu vacancies can be considered as 
the removal of CuO from the diaboleite structure. The 
local transformation is written as follows: 

Pb 2 [Cu(OH) 4 ]Cl 2 — -> [Pb 2 (OH) 2 ]Cl 2 (H 2 0) + CuO. 

This formal equation reflects the fact that the complex 
anion [Cu(OH) 4 ] 2_ , with four OH groups attached to 
the Cu atom, is transformed into the complex cation 
[Pb2(OH)2] 2+ , where the two remaining OH groups are 
attached to Pb. To evaluate the energetics of Cu vacan- 
cies, we use a similar equation and replace the fictitious 
(and presumably unstable) [Pb 2 (OH) 2 ]Cl 2 (H 2 0) com- 
pound with the aforementioned metastable diaboleite 
having 25 % of Cu vacancies (Fig. [3| . Using fully relaxed 
atomic configurations of this Cu-deficient diaboleite, pris- 
tine stoichiometric diaboleite, and CuO (tenorite)pS we 
obtain the energy of +0.27 eV/f.u. (about 10 kJ/mol) 
for the formation of 25 % of Cu vacancies. Although 
rather large, this additional energy can be partially tol- 
erated by the entropy term and by chemical potentials 
of different ions under specific growth conditions. Note 
that Pb2Cu(OH)4Cl2 is prepared in water solution^ and 
likely follows a similar route of crystal growth in the nat- 
ural environment. 

Finally, we consider the problem of stacking disor- 
der. The structure of diaboleite is formed by [Pb2Cl2] 
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slabs that are linked to the [Cu(OH)4] units via the 
Pb-0 bonds. On the opposite side of the slab, the 
connection to the next [Cu(OH)4] magnetic layer is re- 
stricted to weak H- • • CI hydrogen bonds (Fig. |2J top 
left panel). Each CI atom is bonded to four hydrogen 
atoms of the neighboring unit, with the H-Cl distances 
as large as 2.37 A. This value approaches the upper limit 
of O-H- • • CI distances reported in the literature!^ The 
weak interlayer interaction enables regular shifts of the 
neighboring [Pb2Cl2][Cu(OH)4] layers, and facilitates the 
stacking disorder. 

In the tetragonal diaboleite structure, the layers are 
stacked on top of each other (displacement vector k = 0) . 
An alternative stacking sequence could be based on 
k = (||0) because the [PD2CI2] slabs and even oxy- 
gen atoms are invariant to this translation, which only 
changes the order of Cu atoms and vacancies (note the 
Pb2CuDX 6 formula discussed in Sec.[l| along with hydro- 
gen atoms. The resulting structure is monoclinic, with 
the twice larger unit cell and the Cm symmetry (Fig. [4]). 
In GGA+/7, it lies only 40 meV/f.u. above the ideal 
tetragonal structure. This small energy difference is com- 
parable to the entropy term and allows for the formation 
of stacking faults that are indeed observed experimen- 
tally. 

To understand the effect of these stacking faults on 
the reflection broadening, we note that the stacking se- 
quence with k = (§|0) leads to a monoclinic unit cell 
with the lattice vectors a' = a + b, b' = a — b, and 
c' = _|_ Cj where a, b, and c are lattice vectors of the 
parent tetragonal unit cell. The a'b' plane of the mon- 
oclinic cell is fixed as the plane of the diaboleite layer. 
However, the d axis can have two different orientations, 
± c. Its position also depends on the precise value 
of the monoclinic angle between a' and d . Therefore, in 
the reciprocal space the direction of the (c')* axis is fixed 
[(c')* _L a',b', i.e., (c')*||c], whereas the direction of the 
(a')* axis is variable. This explains why the (00/) reflec- 
tions remain narrow, although the reflections with large 
h and k substantially broaden. Indeed, the refinement of 
the anisotropic strain broadening yields the largest hkO 
component, which is related to the variable position of 
the (a') 



axisP 



Although our detailed study of the crystallographic is- 
sues is largely motivated by effects observed in the nat- 
ural sample of diaboleite, similar features are likely rele- 
vant to the Pb2Cu(OH) 4 Ci2 phase in general. The lay- 
ered nature of the crystal structure facilitates the stack- 
ing disorder, whereas the formation of Cu vacancies can 
be instrumental in tuning low-temperature magnetism of 
this material. In the following, we report experimental 
data on thermodynamic properties of diaboleite and ana- 
lyze them from the microscopic viewpoint by considering 
the parent tetragonal crystal structure of stoichiometric 
diaboleite along with possible deviations from this ideal 
atomic arrangement. 




-2 
Energy (eV) 

FIG. 5. (Color online) LDA density of states (DOS) for the 
stoichiometric diaboleite. The Fermi level is at zero energy. 



C. Microscopic magnetic model 

The LDA energy spectrum of diaboleite is typical for 
a cuprate compound. The valence band is dominated 
by oxygen 2p states, with sizable contributions of CI 3p 
around —3 eV, and Cu 3d above —2 eV (Fig. [5]). The 
large underestimate of strong electronic correlations leads 
to the spurious metallic energy spectrum in LDA. The 
LSDA+f/ and GGA+U calculations yield the band gap 
of about 2.7 eV in accord with the blue color of diaboleite 
crystals (Fig. [I] middle panel). 

The Fermi level is crossed by a single band (Fig. IfjJ), 
which is of Cu d x 2_ y 2 origin according to the conven- 
tional crystal-field levels of a strongly elongated octahe- 
dron (the in-plane Cu-0 distances are 1.97 A, whereas 
the out-of-plane Cu-Cl distances are 2.55 A and 2.95 A). 
The Cu d x 2_ y 2 band is well reproduced by a tight-binding 
model with the leading nearest-neighbor intraplane hop- 
ping ti = 78 meV and several long-range interactions 
that are all weak (Table [I]). Note that the coupling be- 
tween the magnetic layers is provided by both t± (along 
the c direction) and t' ± (along [101]). Other hoppings are 
below 2 meV and can be safely neglected in the minimal 
microscopic modelP^ 

Introducing the LDA hoppings into a Hubbard model 
with the effective on-site Coulomb repulsion U c g = 
4.5 eVPwe identify Pb 2 Cu(OH) 4 Cl2 as a magnetic in- 
sulator with t/U e g < 0.02. Therefore, perturbation the- 
ory in ti/Ucs enables the evaluation of AFM couplings 
as J AFM = At 2 jU c ft to the lowest (second) order. We 
find J AFM = 63 K and a very weak frustration by the 
second-neighbor coupling (J 2 AFM /J 1 AFM < 0.01). Weak 
interlayer couplings J± and render diaboleite a good 
quasi-2D magnetic system. 

The above estimates are based on the model approach 
that yields the AFM couplings only. However, magnetic 
couplings - especially those that are short-range - may 
also include a FM component. This FM component can 
be taken into account by total-energy calculations (su- 
percell approach) where energies of collinear spin states 
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FIG. 6. (Color online) LDA band structure of stoichiomet- 
ric tetragonal diaboleite (thin light lines) and the fit with 
the tight-binding model for the Cu d x 2_ y 2 band (thick dark 
line). The k path is defined as follows: T(0, 0, 0), X(0.5, 0, 0), 
Af (0.5, 0.5,0), 2(0,0,0.5), #(0.5,0.5,0.5), .4(0.5, 0,0.5) in 
units of 2ir/a, 2n/b, and 2tt/c, respectively. 



TABLE I. Cu-Cu distances (in A), transfer integrals t; 
(in meV), and exchange couplings Ji (in K) in the tetrag- 
onal stoichiometric diaboleite. The AFM components are ob- 
tained as J, AFM = 4tf/E/eff with U c a = 4.5 eV, whereas the 
Ji values are calculated using the supercell approach. The 
superexchange pathways are shown in Figures [T]and[7| 
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are mapped onto the Heisenberg model to yield full ex- 
change couplings Ji (Tablell|. Then the FM component is 
evaluated as jf M = Ji — J s . In the case of diaboleite, 
this procedure only slightly changes the microscopic sce- 
nario. According to Table [IJ J\ remains the leading in- 
teraction in diaboleite, whereas other couplings are weak 
and AFM. 

To elucidate the origin of magnetic couplings in dia- 
boleite, we calculated the Wannier functions (WF) based 
on the Cu d x i_ y i orbital character. Each WF (Fig. [7]) 
includes the Cu d x 2_ y 2 orbital along with the crp or- 
bitals of oxygen atoms. The superexchange pathway 
of Ji involves the 90° overlap of the neighboring WFs 
(Fig. [7]) , which generally leads to a FM contribution to 
the exchange.^ Indeed, we find a sizable FM compo- 
nent J™ = Ji - J X AFM = -25 K (Table 0. The second- 
neighbor coupling J 2 features a more favorable 180° path- 
way. However, the 0-0 distance of 4.43 A is too large 
for an efficient overlap of the oxygen orbitals, and the 
resulting coupling is very weak. 

Finally, we evaluate magnetic couplings for the stoi- 
chiometric diaboleite with the displaced layers (Fig. [4]) 



FIG. 7. (Color online) Cu d x 2_ y 2 -based Wannier functions, 
and the scheme of Cu-O- • • O-Cu superexchange pathways in 
diaboleite. Different colors show different signs of the Wannier 
function. Larger (green) and smaller (orange) sheres denote 
the Cu and O atoms, respectively. The nearest-neighbor cou- 
pling Ji features the 90° orbital overlap and shorter O • • O 
distances (3.13 A), whereas J2 corresponds to the 180° over- 
lap and longer O- • ■ O distance (4.43 A). 

and for the diaboleite with \ of Cu vacancies (Fig. [3J). 
The formation of vacancies does not disrupt the couplings 
between the remaining Cu sites: the intralayer hoppings 
t\ are in the range 74 — 79 meV, which is very close to 
t\ = 78 meV for the stoichiometric compound. The effect 
of stacking faults is also weak. However, the change in the 
stacking sequence should modify the regime of interlayer 
exchange by replacing Jj_ and J^ with a single coupling 
J", which is somewhat weaker than the interlayer cou- 
plings in the tetragonal structure (compare ~ 3 meV 
with t±_ and t' ± of ±10 meV in Table The role of in- 
terlayer couplings in diaboleite is further discussed in the 
next section. 



D. Thermodynamic properties 

Magnetic susceptibility of diaboleite shows a maximum 
at T max ~ 32 K followed by a kink at T N = 11 - 12 K 
(Fig. [8]). This kink is somewhat field-dependent (Fig. [9]) 
and indicates the onset of long-range magnetic order. 
Above 200 K, the susceptibility follows the Curie- Weiss 
law 

C 



with xo = — 2 x 10~ 4 emu/mol, C = 0.512 emu K/mol, 
and = 35 K. The sizable temperature-independent dia- 
magnetic contribution is likely related to the contribution 
of the sample holder and core diamagnetism. The C value 
yields the effective magnetic moment // c ff = 2.02 /j,b, 
which is notably higher than 1.73 hb expected for spin-|. 
The deviation of n e g from the spin-only value implies 
g ~ 2.32, which is rather high but still in the reason- 
able range for Cu +2 compounds (see, e.g., Ref. 49). The 
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FIG. 8. (Color online) Magnetic susceptibility of diaboleite 
measured in the applied field of 5 T (circles) and QMC fits 
of three spin models (lines): the stoichiometric square lattice 
and the square lattices with 5% and 10% of vacancies. The 
arrow denotes the Neel temperature Tm — 11.5 K at 5 T. 
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FIG. 9. (Color online) Field dependence of the Neel temper- 
ature in diaboleite, as seen from the magnetic susceptibility 
and specific heat (inset) data. Lines are guide for the eye. 
Arrows denote the transition temperatures (see text for de- 
tails). 



positive 9 value implies predominantly AFM exchange 
couplings. 

The susceptibility maximum can be approximated by a 
QMC fit for the purely 2D Heisenberg model of the spin- \ 
square lattice (solid line in Fig. [8]). We find J = 35 K, 
g = 2.34, and xo = — 2 x 10~ 4 emu/ mol in excellent agree- 
ment with the high-temperature Curie- Weiss fit. How- 
ever, the square-lattice model fails to describe the data 
below the susceptibility maximum. In this region, the 
fit can be improved by considering a diluted square lat- 
tice with randomly distributed vacancies. The dilution 
of spin-i sites with vacancies increases the susceptibility 
below T ,max and eventually blurs the maximum. We find 
the best fit for the lattice with 5 % of vacancies, J = 37 K, 
and nearly unchanged g and xo parameters. This result 
perfectly matches J x = 38 K from DFT (Table The 
concentration of vacancies is well in line with the XRD 
data that yield about 6% Cu deficiency (Sec. |III A| . 

The shape of the susceptibility curve below T max is 
a strong evidence for the formation of vacancies in the 
magnetic layers of diaboleite. Although 3D correlations 
(interplane couplings) might have similar effect on the 
susceptibility, their contribution is not strong enough to 
explain the experimental data. For example, weakly cou- 
pled square planes with an effective interlayer coupling 
J±*/Ji — 6.5 x 1CP 3 , which reproduces the magnetic 
transition temperature in diaboleite (see Fig. 10), show 



nearly the same susceptibility as the single square plane. 

The antiferromagnetic ordering in diaboleite manifests 
itself by a kink of the magnetic susceptibility (Fig. [9]). 
The Neel temperature increases with the applied field. 
The magnetic transition is not well seen in zero-field spe- 
cific heat, although a close examination of the C p /T data 
at 11 — 12 K spots a tiny feature at 10.7 K (see the inset of 
Fig. 10). Otherwise, the temperature dependence of the 




zero-field specific heat is smooth up to at least 100 KPS 
A magnetic field enhances the transition anomaly that 
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FIG. 10. (Color online) Field-temperature phase diagram de- 
rived from QMC simulations for weakly coupled square planes 
with an effective interlayer coupling J± /Ji = 6.5 x 10~ 3 . Ex- 
perimental values of TW are shown according to the magnetic 
susceptibility (x) and specific heat (C p ) data. Left (colored) 
and right (white) parts are the magnetically ordered and para- 
magnetic phases, respectively. The inset shows the specific 
heat measured in zero field. 



becomes well visible at 6 T and 9 T (see the inset of 
Fig.|9j). 

The broadening of the transition anomalies, which is 
inevitable in powder samples, prevents us from the pre- 
cise evaluation of Tjy in different magnetic fields. Nev- 
ertheless, the data in Fig. [9] show clearly that Tjv in- 
creases with the applied field, and this trend holds up to 
at least 9 T. Using middle points of the local minima of 
the susceptibility, as well as the kinks in the specific heat 
data, we tentatively evaluate the field dependence of Tjv 
(Fig. 10). The steady increase in Tjv is very typical for 
spin- j quantum magnets and resembles the field-induced 
behavior observed in other square-lattice systems with 
spin. iJEEZI 



To analyze the field dependence of Tjv, we first con- 
sider the magnetic transition in zero field. The 3D AFM 
order is driven by the interlayer couplings J± and J',. 
These couplings are comparable in strength (see Tabled) 
and frustrated, because they favor different interlayer 
spin configurations: G-type (antiparallcl spins on nearest 
neighbors along c) and C-type (parallel spins on near- 
est neighbors along c), respectively. However, the DFT 
evaluation of J± and J' ± comes with a caveat that the 
couplings on the order of 1 K are very difficult to esti- 
mate precisely, and a careful experimental work, such as 
inelastic neutron scattering at very low energies, would 
be necessary to get accurate values of J± and J' ± . Here, 
we argue that the experimental data are in agreement 
with an effective non-frustrated interlayer exchange that 
might be a combination of J± and J' ± . 

We have studied the AFM ordering in weakly coupled 
square planes by QMC simulations of a 3D spin model 
with the uniform interlayer coupling J^ ff (see Refs. l37l 
and l50l for similar studies). The J^ ff parameter was ad- 
justed as to match T/v ~ 10.7 K in zero field. The result- 
ing estimate of Jf/Ji = 6.5 x 1(T 3 (i.e., ~ 0.2 K) 
is somewhat smaller than Jj_ and J' ± (TableTTl) and may 
reflect their partial compensation. Using this value of 
J^ ff , we are able to reproduce both T/v in zero field and 
the field dependence of T/v (Fig. 10). This way, the ex- 



perimental data available so far are in agreement with 
the spin model of weakly coupled square planes without 
any substantial frustration. Regarding the role of Cu va- 
cancies and stacking faults, their effect on the interlayer 
exchange is likely minor and can not be distinguished 
using the experimental data available so far. 

The physical mechanism behind the increase in T/v is 
the suppression of quantum fluctuations by the applied 
magnetic field. In higher fields, this effect is countered 
by the increasing tendency toward parallel spin arrange- 
ment (see, e.g., Ref . l37|) . However, the AFM couplings in 
diaboleite are strong enough to compete with the field, 
so that the T/v increases linearly up to at least 9 T. 



IV. DISCUSSION AND SUMMARY 

Our study of diaboleite raises several problems that 
may have implications for other quantum magnets and 
for theoretical models of quantum magnetism. First, 
our anticipation of the sizable second-neighbor coupling 
J2 was not confirmed in the experimental and compu- 
tational study. While diaboleite indeed bears struc- 
tural similarity to PbV03, which is the strongly frus- 
trated square- lattice compound with sizable J2> de- 
tails of their electronic structures are notably different. 
In PbVC>3, vanadium polyhedra are directly connected 
to each other, thus leading to relatively short pathways 
for Ji (V-V distance of 3.88 A) and J 2 (5.49 A). The 
diaboleite structure features isolated CuO-4 plaquettes, 
hence both Ji and J 2 correspond to much longer Cu-Cu 
distances of 5.88 A and 8.32 A, respectively. 



The substantial elongation of the metal-metal dis- 
tances in diaboleite could be mitigated by a different 
symmetry of the magnetic 3<i orbital. In contrast to 
the half- filled d xy orbital of V +4 , the d x 2_ y 2 orbital of 
Cu +2 features the strong er-overlap with the oxygen p or- 
bitals. This way, the Wannier functions are extended by 
large oxygen contributions ( "tails" ) that can overlap and 
eventually lead to strong couplings even for very long su- 
perexchange pathways. In diaboleite, this mechanism is 
indeed operative for J\ , where the 0-0 distance of 3.13 A 
results in a sizable coupling of 35 K despite the 90° over- 
lap, which is generally deemed unfavorable for the AFM 
coupling. Regarding J 2 , the 0-0 distance of 4.43 A ap- 
pears to be too long for any appreciable magnetic inter- 
action, even though the orbitals form the favorable 180° 
geometry (Fig. [7| . 

The second interesting aspect of the diaboleite mag- 
netism pertains to experimental signatures of the mag- 
netic transition. While the susceptibility measurements 
unequivocally demonstrate the 3D magnetic ordering be- 
low T/v ~ 11 K (Fig. [9j), heat capacity measured in zero 
magnetic field shows only a tiny feature that, when exam- 
ined on its own, may be well considered an artifact (see 
the inset of Fig. 10 ). The lack of a conspicuous transition 



anomaly is caused by the small amount of the magnetic 
entropy available at T/v. Using QMC simulations, we es- 
timate that in zero field only 7 % of the magnetic entropy 
can be released at T/v . The main release of the entropy is 
observed at T ~ J as a broad maximum in the temper- 
ature dependence of the magnetic specific heatP^^l An 
external magnetic field shifts the entropy from the broad 
maximum at T ~ J toward the transition anomaly at 
T/yP^ This way, the transition anomaly becomes visible 
and can be tracked with heat-capacity measurements. 

Similar effects of low magnetic entropy at T/v and its 
dramatic increase in the applied field have been reported 
for the square-lattice magnets with weak exchange cou- 
plings and relatively low saturation fieldsP^ESl Although 
diaboleite features the stronger exchange and the high 
saturation field (expected at H s = AJigns — 120 T), 
even the low fields of 6 - 9 T (H/H s = 0.05 - 0.075) are 
sufficient to render the transition anomaly visible. Our 
experimental results are in line with earlier reports on 
other spin-| materials. For example, Lancaster et a/P^l 
claimed that quasi-2D systems with weak interlayer ex- 
change do not show transition anomalies in the spe- 
cific heat, whereas Sengupta et al^ demonstrated the 
suppression of the transition anomaly depending on the 
size of the interlayer couplings. Our data emphasize 
the importance of heat-capacity measurements in high 
magnetic fields, because an external field enhances the 
anomaly and facilitates an experimental observation of 
the magnetic transitions. This result may be relevant 
to t he co ntroversial magnetic transition in PbV03 as 
well™ 

Diaboleite is a rare material showing magnetism of the 
diluted square lattice with spin- 4. We have proved the 
formation of Cu vacancies by XRD and chemical analy- 
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sis, and observed the dilution effect in magnetic suscep- 
tibility measurements (Fig. [8]). If the concentration of 
Cu vacancies could be varied (e.g., by crystal growth in 
Cu-deficient environment), diaboleite becomes an inter- 
esting system with a remarkably "clean" mechanism of 
dilution. Presently, most studies of the diluted square 
lattice focus on Zn- and Mg-doped La2Cu04, where the 
doping with non-magnetic atoms leads to a frustration of 
the spin lattice by the second-neighbor coupling J 2 ! 53 l 54 l 
In diaboleite, this effect is suppressed, because the Cu- 
Cu distances are very large and prevent any long-range 
superexchange interactions (Fig. Y7n . The diluted spin- \ 
systems provide a unique opportunity to study the inter- 
play of disorder and quantum fluctuations EUE3 there- 
fore, further studies of diaboleite may be insightful. 

More generally, an in-depth physics research performed 
on natural samples is a delicate compromise between the 
potential complexity of the sample and the tantalizing 
opportunity to find unexpected structural and physical 
effects. Natural samples are prone to contamination by 
foreign phases and elements. A meticulous sample char- 
acterization is an essential part of the work that often 
ends when the sample is found unsuitable, owing to a 
contamination by magnetic impurities. However, sam- 
ple imperfections can be advantageous as well. Some- 
times they disclose interesting effects that are neither ex- 
pected nor observed in synthetic samples. We hope that 
the study of Cu-based minerals will further contribute to 
the experimental and theoretical work on quantum mag- 
netism. Detailed characterization of the chemical compo- 
sition, crystal structure, and physical properties should 
also advance our understanding of intricate natural pro- 
cesses that lead to the formation of these interesting and 



aesthetically beautiful materials. 

In summary, we have studied the crystal structure and 
magnetism of the diaboleite mineral. The formation of 
Cu vacancies, which are tolerated by the crystal struc- 
ture, leads to the dilution of the spin lattice. This effect 
is probed directly via chemical analysis and XRD, as well 
as indirectly via magnetic susceptibility measurements 
that show clear signatures of the dilution. The concen- 
tration of vacancies is about 5%. Consistently, thermo- 
dynamic properties of diaboleite are well described by the 
spin model of the diluted square lattice with the nearest- 
neighbor exchange J\ ~ 37 K and 5% of Cu vacancies. 
Weak interlayer couplings trigger the long-range mag- 
netic order below Tn — 10.7 K in zero field. 
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Supplementary information for 
Square-lattice magnetism of diaboleite Pb 2 Cu(OH) 4 Cl 2 
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FIG. 11. Rietveld refinement of the room-temperature high-resolution XRD data. Ticks 
show the reflection positions, from bottom to top: diaboleite (D), phosgenite (P), and 
wherryite (W). Note that for the diaboleite phase we had to refine the anisotropic strain 
broadening parameters S hM (in 1CT 2 deg 2 A" 4 ): S 4ao = 5 004 = -0.033(2), S 220 = 0.119(3), 
S202 = S022 = 0.162(3), S004 = 0. The definition of Shki is according to the Jana2006 
program, see also [P. W. Stephens, J. Appl. Cryst. 32, 281 (1999)]. 



TABLE II. Atomic positions, atomic displacements parameters Ui ao (in 10 -2 A 2 ), and occupation factors g 
refined from the powder data on the natural sample of diaboleite. Although the accuracy of the structural 
parameters is probably lower than in [M. A. Cooper and F. C. Hawthorne, Canad. Mineralog. 33, 1125 
(1995)], where single-crystal data were used, our refinement provides information on peculiarities of the 
particular sample used in the present work. Error bars are obtained from the Rietveld refinement. 



Atom 


Wyckoff position 


x/a 


v/b 


z/c 




9 


Pb 


2c 





1 

2 


0.7298 a 


1.2(1) 


1.0 


Cu 


1« 








0.008(3) 


1.0* 


0.942(4) 


Cll 


la 








0.464(3) 


0.1(l) c 


1.0 
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0.1(l) c 
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1.3(1) 
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Id 
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a The z coordinate of the Pb atom is fixed in order to keep the position of origin 
6 Only the occupancy factor was refined 

c The atomic displacements parameters of Cll and C12 were refined as a single parameter 
d The position of hydrogen was not refined 




FIG. 12. Zero-field heat capacity of diaboleite. 



